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We present the current controlled motion of dynamic soliton embedded in spin wave background 
in ferromagnetic nanowire. With the stronger breather character we get the novel magnetic rogue 
wave and clarify its formation mechanism. The generation of magnetic rogue wave is mainly arose 
from the accumulation of energy and magnons toward to its central part. We also observe that 
the spin-polarized current can control the exchange rate of magnons between envelope soliton and 
background, and the critical current condition is obtained analytically. Even more interesting is 
that the spin-transfer torque plays the completely opposite role for the cases of below and above 
the critical value. 
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I. INTRODUCTION 

In ferromagnet nanowires, the deviation of magnetiza- 
tion from the ground state results in the excitation of 
spin waves. Their attractive interaction and instabilities 
contribute to the existence of topological and dynamic 
solitons. The dynamics of magnetization in a soliton can 
be well described by the famous Landau-Lifshitz-Gilbert 
equation [1]. 

The topological soliton, now namely domain wall 
(DW), connecting two vacua describes an inhomogeneous 
state of magnetization and it cannot be reduced to the 
ground state by any finite deformation Q. So it has 
technological application in race track memory The 
field-driven DW motion has been studied extensively in 
thin films [4] and recently in magnetic nanowires (EHIoj. 
These studies show that the DW travelling speed in mag- 
netic materials conventionally depends on the strength 
of a constant magnetic field as the external filed is below 
the Walker-breakdown value [ll|, [l2| . Above this critical 
value, the DW develops a variational complex internal 
structure which results in temporal oscillation of DW 
velocity [!, [I, 0, [l| . Recently, many studies demon- 
strate that the spin-polarized current can cause many 
unique phenomena for magnetization motion, which at- 
tributes to spin-transfer effect [15[. With this considera- 
tion the modified Landau-Lifshitz-Gilbert equation p~Gl — 
fl8| including spin-transfer torque has been derived to 
describe such current-induced magnetization dynamics. 
With the remarkable experimental successes measuring 
the motion of a DW, considerable progresses have been 
made to understand the current-induced DW motion in 
magnetic nanowires [l7l-[23|. This current induced mag- 
netization motion possesses of different character from 
that driven by an magnetic field, and the spin-transfer 



torque can induce not only precession but also damping 
role for the motion of magnetization. 

On the other hand, the dynamic soliton describes the 
localized states of magnetization which can be reduced 
to a uniform magnetization by continuous deformation, 
so that the excited ferromagnet makes a transition to 
the ground state. Therefore, a dynamic soliton is some- 
times said to be topologically equivalent to the ground 
state. Also, the motion of dynamic soliton is of topic 
research in confined ferromagnetic materials (23-l27j. es- 
pecially with the generation and detection of magnons 
excitation [28| in a magnetic multilayer. Driven by the 
adiabatic spin-transfer torque, the dynamic soliton solu- 
tions for isotropic case [29| and uniaxial anisotropic case 
[30| are investigated carefully, where the corresponding 
solutions show the characteristic breather behavior for 
magnetization motion in ferromagnetic nanowire. How- 
ever, the current driven motion of the dynamic soliton 
are not well explored carefully. 

In this paper, we present the generation of a novel 
magnetic rouge wave solution under stronger breather 
characters in a perpendicular anisotropic ferromagnetic 
nanowire driven by spin-transfer torque, which is similar 
to optical rouge wave in fiber [3l|, |32|. The formation 
mechanism of magnetic rouge wave is clarified carefully. 
In this process, we observe that the spin-polarized cur- 
rent can control the magnons exchange rate between the 
envelope soliton and background, and the critical current 
condition is obtained analytically. Moreover, the spin- 
transfer torque play the completely opposite role for the 
cases of below and above the critical value. 



II. CURRENT DRIVEN INTERACTION OF 
SPIN WAVE AND A DYNAMIC SOLITON IN 
FERROMAGNETIC NANOWIRE 



* Corresponding author, E-mail: zdli2003@yahoo.com 



When a spin-polarized electric current flows the ferro- 
magnetic nanowire, the localized magnetization dynam- 



ics can be described by the modified Landau-Lifshitz 
equation [l?], [H| with the spin- transfer torque, 
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where 7 is the gyromagnetic ratio, a is the Gilbert damp- 
ing parameter, and H e ff represents the effective magnetic 
field including the exchange field, the anisotropy field and 
the external field. For a perpendicular anisotropic ferro- 
magnetic nanowire, the effective magnetic field takes the 
form H eff = (2F/M 2 ) d 2 M/dx 2 + [(H K /M S - 4tt) M z + 
H ex t]e z , where is the exchange constant, Hk denotes the 
energetic anisotropy coefficient, H ext is the applied exter- 
nal magnetic field, and e z is the unit vector along the z- 
axis. The last term r s in Eq. (pQ) denotes the spin- transfer 
torque, which describes the injected current can be po- 
larized and produces a torque acting on the local mag- 
netization. This adiabatic spin-transfer torque is of the 
form r s = bj(dM/dx) [Hj], where bj = Pj e ^ B / (eM s ), 
P is the spin polarization of the current, j e is the electric 
current density and flows along the x direction, \ib is the 
Bohr magneton, e is the magnitude of electron charge. 

Introducing the normalized magnetization, i.e., m = 
M/Ms, we can simplify Eq. (pp) as the dimensionless 
form 

dm , dm . dm . , 

— = -mxh eff + amx— + Aj — , (2) 

where h e ff = H e ff/M s and Aj = bjto/lo. The time t 
and space coordinate x have been rescaled by the char- 
acteristic time to — 1/ [7 (Hk — 47rM s )] and length lo = 
\/2Fj [(Hk — 4:7rM s ) M s ], respectively. As discussed in 
our pervious work [30|, Eq. (|2j) admits the excited states, 
i.e., spin wave and a dynamic soliton. For the perpendic- 
ular anisotropic ferromagnetic nanowire, these two types 
of excited state denote the small deviation of magnetiza- 
tion from the ground state. Therefore, a reasonable com- 
plex function q can be introduced to replace the compo- 
nents of normalized magnetization [2|, i.e., q = m x + im y 
and m 2 = 1 — \q\ 2 . Under the long- wavelength approx- 
imation Q and without damping, Eq. (j2j) becomes the 
integrable nonlinear Schrodinger equation 



.dq 



dx 2 



1 I I 2 



■iAj 



dx 



(3) 



where u = 1 + H ext / (H k - 4ttM s ). 

It is easy to find two basic solutions of Eq. (|3]). 
One is q = 0, which corresponds to the ground state 
m = (0,0,1). The other solution is spin wave, i.e., 
q = A c e~ l<K - kcX ~ UJct ^ ) with uj c and k c being the dimension- 
less frequency and wave number, respectively. In per- 
pendicular anisotropic ferromagnetic nanowire, the inter- 
action of spin waves is attractive and their instabilities 
lead to macroscopic phenomena, i.e., the appearance of 
a spatially localized magnetic excited state (topological 
or dynamic soliton). As shown later, the interaction of 
spin wave and a dynamic soliton can result in the novel 
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FIG. 1: (Color online) Breather velocity V q and width W vs 
the amplitude A c and wave number k c of spin wave. The 
other parameters are ui = 0.2, isi = —0.3 and Aj = 0.2. (a) 
and (d): The soliton velocity V q and width W vs with A c and 
k c . The red star signs the location of the rogue wave, (b) and 
(c) The dependence of velocity on A c and k c . (e) and (f) The 
dependence of velocity on A c and k c . 



magnetic rogue wave which can be realized by adjusting 
the relation of wave number and amplitude for spin wave 
and a dynamic soliton. To this purpose, one should get 
the solution on the spin wave background. 

By employing Darboux transformation [30|, [33| and 
performing a tedious calculation, we obtain such solution 



q = e i *[A c + 2u 1 (Ai+iA 2 )/A], 



(4) 



where tp = u c t—k c x, A = cosh 6+a cos /3, Ai = a cosh 6+ 
cos /?, A2 = 6sinh# + csin /?, here = DRX J r(D5)Rt J rXQ 1 
P = D T x + (DS)jt- t , a = 2A C L R /(\L\ 2 + A 2 ), b = 
2A C L T /(\L\ 2 + A 2 ), c=(\L\ 2 - A 2 C )/(\L\ 2 + A 2 C ), and the 
subscript R and / represent the real and imaginary part, 
respectively. The other parameters are L = —ik c — D — \, 
D = [(ik c + A) 2 - A 2 ] 1 / 2 , and S = -iX - k c + Aj, where 
A = ui + iv\. Here xo, to, A c , fc c , u\, and v\ are the real 
constants, and without loss of generality we assume that 
A c and u\ are non-negative real constants. 

The solution in Eq. (jlj) exhibits a dynamic soliton so- 
lution embedded in spin wave background. Two types of 
excited states (spin wave and a dynamic soliton) can be 
recovered from Eq. (j4]): a) As the spin wave amplitude 
and wave number vanishes, namely A c = k c = 0, we 
can get a dynamic soliton, q\ — 2uie l ^ 1 ^ UJct " ) / cosh #1, 
with 0i = u\x + u\ (2v\ + Aj) t + xq and f3\ = v\x — 
(u 2 — v\ — Ajv\) t — to- This solution corresponds to 
the dynamic precession of magnetic soliton on the ground 
state background, where the soliton moves with the ve- 
locity v = —2v\ — Aj, and the components m x and m y 
precess around m z with the amplitude A q = 2u\j cosh#i 



and the frequency uj q 



Ajv x 



It clearly 



demonstrates that the spin current term Aj can change 
the dynamic soliton velocity and the precessional fre- 
quency. In addition, the solution qi represents in fact the 
static magnetic soliton with three integrals of the motion 
and the uniform magnon density \qi\ = 2u\ along the 
direction of soliton propagation, b) When the amplitude 
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FIG. 2: (Color online) The evolution of magnetic soliton in 
the limit processes u± — >• A~ (a) and i^i — >• A J (b). Other 
parameters are A c — 0.2, Aj — 0.2, fc c = —vi — 0.003, and 
xo = to = 0. 



of dynamic soliton vanishes, namely u± = 0, Eq. (j4j re- 
duces to spin wave solution q = A c e~ l<K - kcX ~ UJct \ in which 
the magnon density is also constant. 

From Eq. ((4]), we observe that the solution commonly 
exhibits a breather character and a time periodic mod- 
ulation for soliton amplitude, which in fact denotes the 
interaction between the localized process of spin wave 
background and the periodization process of magnetic 
soliton. The properties of the soliton solution are char- 
acterized by the slope direction Vq = Drx + {D5)nt, the 
breather propagation velocity v q = —(D5)r/Dji and the 
soliton width 1/Dr. With the expressions of D and 
we find that the soliton velocity and width are modu- 
lated by the amplitude A c and wave number k c of spin 
wave as shown in Fig. 1(a) and 1(d). The absolute value 
of soliton velocity and soliton width increases with in- 
creasing A c as shown in Fig. 1(c) and 1(f). Moreover, 
as the modulation parameter A c increases continuously, 
the soliton velocity will be affected mainly by the value 
of k c near —v\. When k c = —v\, the soliton velocity and 
width is maximal, respectively, as shown in Fig. 1(b) 
and 1(e). This phenomena in fact ascribes to magnonic 
spin-transfer torque [34[ which denotes the transfer of 
spin angular momentum from spin wave background to 
a dynamic soliton, and it will be discussed elsewhere in 
detail. On the basis of Eq. (jlj, as the modulation pa- 
rameter u\ approaches A C: some novel properties will be 
presented. Magnetic rouge wave solution will be excited 
in ferromagnetic nanowire which leads to some fantastic 
phenomena. 



III. NOVEL MAGNETIC ROGUE WAVE 

In terms of our previous discussion [30j for the solu- 
tion in Eq. (j4]), we have known that the critical point 
u\ = A c forms a dividing line between the modulation 
instability process (u± > A c ) and the periodization pro- 
cess (u\ < A c ) under the condition v\ = — k c . It leads 




FIG. 3: (Color online) The magnon density distribution 
against the background for the different parameter u\ , which 
ranges from 0.09 to 0.29 in 0.05 steps, (a) The magnon density 
distribution against the background for the excited formation 
of magnetic rogue wave. Other parameters are A c = 0.2, 
A j = k c = 0.1 and xq = to = 0. 



to the different physical behavior how the breather char- 
acter depends strongly on the modulation parameter u\ 
as shown in Fig. 2. Fig. 2 shows two different asymp- 
totic behavior in the limit processes u\ —> A~ (and A£) 
under the condition v\ = — fc c , respectively. The for- 
mer demonstrates a spatial periodic process of a soliton 
and we observe that the spatial separation of adjacent 
peak and each peak value increases rapidly as the mod- 
ulation parameter u\ approaching A c , respectively. The 
latter shows a localized process of the spin-wave back- 
ground along the slope direction K$ = —(D5)i/Dj for 
to = 0. In this case, the temporal separation of adjacent 
peak also increases rapidly as the modulation parame- 
ter u\ approaching A C1 while each peak value takes the 
rapid decrease. Even more interesting, in the limit case 
of u\ —> A c , i.e., the pentagram sign indicated in Fig 1(a) 
and (d), we get the novel magnetic rogue wave 



4(l-i*Ag) 
t 2 A 2 c r] + 2txA 2 c ( + e 



(5) 



where 77 = A 2 j + A 2 C + Ak 2 - AAjk c , ( = Aj - 2k c , and 
e = 1 + x 2 A 2 . Eq. (j5]) shows the typical rogue wave 
feature that the magnons accumulated from spin wave 
background converge a single hump with the critical am- 
plitude Aq = 3A C . It implies that the localization wave 
is captured completely at x = and t = by spin wave 
background. As shown in Fig. 3, the realization of mag- 
netic rogue wave in Eq. (j5j) attributes to that the aggre- 
gation of magnons gradually increases as u\ approaches 
A~ , while decreases gradually with u\ approaches A£. 
The temporal localization magnetic rogue wave is excited 
as Ui — >• A c with the magnon density peak \q\ 2 = 9 A 2 . 

In order to investigate deeply the properties of mag- 
netic rogue wave in ferromagnetic nanowire, the analysis 
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FIG. 4: (Color online) (a)-(e) The formation region in space 
(x, t) for magnetic rogue wave with different current. The 
parameter Aj ranges form to 0.8 in 0.2 steps, (f) The 
nonuniform exchange of magnons between rogue wave and 
background for the different spin current. The inset figure in 
(f) denotes the maximal accumulation (or dissipation) process 
for the critical current value Aj — 2k c . Other parameters are 
A c 0.2 and k c 0.2. 



for the magnon density distribution against the back- 
ground plays a major role, which is defined by the quan- 
tity p q (x,t) = \qi (x,t)\ 2 -\q! (x = ±oo,£)| 2 . In Fig. 3 we 
plot the evolution of magnon density distribution for the 
breather solution in Eq. (|4]), which ultimately demon- 
strates the formative mechanism of magnetic rogue wave 
solution in ferromagnetic nanowire. With the modulation 
parameter u\ approaches A c , the magnons in the back- 
ground gradually gather toward to each individual cen- 
tral part and the envelope becomes sharper. Specially, as 
shown the inset (a), the critical peak will appear under 
the condition u\ —> A c , i.e., the magnetic rogue wave is 
excited. 

With Eq. (j5j) and the quantity p = \Qi(x,t)\ 2 — 
\Qi (x = ±oo,£)| 2 we obtain the magnon density distri- 
bution in a magnetic rogue wave 

ri-r 2 



p = 8Ai 



(r! + r 2 f 



(6) 



where T 1 = 1 + t 2 A A c , T 2 = A 2 c (x + t (Aj - 2k c )f , and 
Eq. (j6j) implies the integral p(x,t)dx = for arbi- 
trary time. From the condition pg (±1/A C: 0) = 0, we 
can define the spatial width of the hump part in rogue 
wave as 2/A c . A detail calculation shows that at a fixed 
time the loss of magnons in background completely trans- 
fer to hump, i.e., the area relation Si + S 2 = S3. These 
results clearly illustrate that the generation of magnetic 
rogue wave with stronger breather character is mainly 
arose from the gathering energy and magnons from the 
background toward to its central part, and the loss of 
magnons in background completely transfer to the hump 
part in magnetic rogue wave. 

The other interesting fundamental problem is that 
how rogue wave gather magnons and energy toward to 
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FIG. 5: (Color online) (a) Evolution of the numerical mag- 
netic rogue wave, (b) Comparison of magnon density between 
the numerical solution (in solid curve) and ideal solution (in 
V). Other parameter are A c = 0.2, u\ — 0.19, xo = 0.72, 
Aj 0.01, k c 0.2 and uj 0.2. 



its central part from the background. This can be ex- 
plained by the quantity 5(x,t) = lim/ Q ^±oo \Qi (x,t) — 
Qi (x = lQ,t)\ 2 . With Eq. © we obtain 



S(x,t) = WA 2 C 



(ri + r 2 ) 



2 ' 



(7) 



which denotes the nonuniform exchange of magnons be- 
tween rogue wave and background for the different spin 
current as shown in Fig. 4. From Eq. ([?)) we find 
the spin current can control the accumulation and dis- 
sipation rate of magnons, and there is a critical current 
condition, i.e., Aj c = 2k c . Below the critical current, 
the magnons exchange decreases with the increasing cur- 
rent term Aj. However, the magnons exchange is ac- 
celerated with the increasing current above the critical 
value. The roles of spin- transfer torque are completely 
opposite for the cases below and above the critical cur- 
rent which is shown in Fig. 4 (f). From Fig. 4(a) to 
4(e) we see that the magnetic rogue wave can be cre- 
ated in the different direction for (x, t) space, which as- 
cribes to the nonuniform exchange of magnons between 
rogue wave and background tuning by the spin-polarized 
current. When Aj = 2fc c , the time of magnons ac- 
cumulation (or dissipation) attains its maximum. As 
shown the inset figure of Fig. 4 (f), i.e., the integral 
£(x,t) = J^S{x,t)dx = 8ttA c /(1 + t 2 A A c y/ 2 , the 
magnons in background accumulate to the central part 
when t < 0. It leads to the generation of a hump with 
two grooves on the background along the space direction 
and the critical peak of the hump can occur at t = 0. In 
contrast, when t > 0, the magnons in the hump start to 
dissipate into the background so that the hump gradually 
decay. The magnetic rogue wave disappears ultimately, 
and it verifies the rogue wave is only one oscillation in 
temporal localization and displays a unstable dynamic 
behavior. 

With some specific initial conditions, the excitation of 
rogue wave can be recovered by means of the numerical 
simulation. In order to understand such process of mag- 
netic rogue wave, we choose the initial value of solution, 



5 



which can be approximated by 

q 0, 0) = + ex cos (p ± ) exp (i<p) , (8) 

where p = [2ki (k>i — iu\) — A 2 ] /A c , e = exp (— xq), x — 
^K\U\ (k\ — iu\) /A 2 C and ipi = K\x with = \j A? c — u\. 
By direct numerical simulations, we find that the solu- 
tion of the initial value problem in Eq. (|3]) with initial 
condition Eq. (j8]) can be well described by the solution 
in Eq. (|5]). Fig. 5(a) is a result of direct numerical sim- 
ulation, which shows that a small periodic perturbation 
with a very small modulation parameter can induce a 
near-ideal rogue wave localization, whose profile is basi- 
cally consistent with the ideal theoretical limit solution 
of Eq. (|5|) as shown in Fig. 5(b). As a result, a small 
initial perturbation with a small modulation can induce 
the generation and breakup of a near-ideal rogue wave. 

IV. CONCLUSIONS 

In summary, we have investigated the formative mech- 
anism of magnetic rogue wave and the properties of 
magnon density in uniaxial anisotropic ferromagnetic 
nanowire driven by spin-transfer torque. Our results 
show that the accumulation of energy and magnons to- 
ward to its central part plays the main role for the gen- 



eration of rogue wave with stronger breathing character 
and rogue wave only appears at the spatial-temporal lo- 
calization. We also display the rouge wave is unsuitably 
and gradually decay because the nonuniform exchange 
rate of energy and magnons which can be controlled by 
the spin-polarized current. A novel critical current condi- 
tion is obtained analytically, and the spin-transfer torque 
plays the completely opposite roles for the case of below 
and above the critical value. 
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